Search for: All records

The millisecond pulsar J1713+0747 underwent a sudden and significant pulse shape change between 2021 April 16 and 17 (MJDs 59320 and 59321). Subsequently, the pulse shape gradually recovered over the course of several months. We report the results of continued multifrequency radio observations of the pulsar made using the Canadian Hydrogen Intensity Mapping Experiment and the 100 m Green Bank Telescope in a 3 yr period encompassing the shape change event, between 2020 February and 2023 February. As of 2023 February, the pulse shape had returned to a state similar to that seen before the event, but with measurable changes remaining. The amplitude of the shape change and the accompanying time-of-arrival residuals display a strong nonmonotonic dependence on radio frequency, demonstrating that the event is neither a glitch (the effects of which should be independent of radio frequency,ν) nor a change in dispersion measure alone (which would produce a delay proportional toν⁻²). However, it does bear some resemblance to the two previous “chromatic timing events” observed in J1713+0747, as well as to a similar event observed in PSR J1643−1224 in 2015.

 

Recently we found compelling evidence for a gravitational-wave background with Hellings and Downs (HD) correlations in our 15 yr data set. These correlations describe gravitational waves as predicted by general relativity, which has two transverse polarization modes. However, more general metric theories of gravity can have additional polarization modes, which produce different interpulsar correlations. In this work, we search the NANOGrav 15 yr data set for evidence of a gravitational-wave background with quadrupolar HD and scalar-transverse (ST) correlations. We find that HD correlations are the best fit to the data and no significant evidence in favor of ST correlations. While Bayes factors show strong evidence for a correlated signal, the data does not strongly prefer either correlation signature, with Bayes factors ∼2 when comparing HD to ST correlations, and ∼1 for HD plus ST correlations to HD correlations alone. However, when modeled alongside HD correlations, the amplitude and spectral index posteriors for ST correlations are uninformative, with the HD process accounting for the vast majority of the total signal. Using the optimal statistic, a frequentist technique that focuses on the pulsar-pair cross-correlations, we find median signal-to-noise ratios of 5.0 for HD and 4.6 for ST correlations when fit for separately, and median signal-to-noise ratios of 3.5 for HD and 3.0 for ST correlations when fit for simultaneously. While the signal-to-noise ratios for each of the correlations are comparable, the estimated amplitude and spectral index for HD are a significantly better fit to the total signal, in agreement with our Bayesian analysis.

 

Abstract Neutron Star Interior Composition Explorer has a comparatively low background rate, but it is highly variable, and its spectrum must be predicted using measurements unaffected by the science target. We describe an empirical, three-parameter model based on observations of seven pointing directions that are void of detectable sources. Two model parameters track different types of background events, while the third is used to predict a low-energy excess tied to observations conducted in sunlight. An examination of 3556 good time intervals (GTIs), averaging 570 s, yields a median rate (0.4–12 keV; 50 detectors) of 0.87 c s −1 , but in 5% (1%) of cases, the rate exceeds 10 (300) c s −1 . Model residuals persist at 20%–30% of the initial rate for the brightest GTIs, implying one or more missing model parameters. Filtering criteria are given to flag GTIs likely to have unsatisfactory background predictions. With such filtering, we estimate a detection limit, 1.20 c s −1 (3 σ , single GTI) at 0.4–12 keV, equivalent to 3.6 × 10 −12 erg cm −2 s −1 for a Crab-like spectrum. The corresponding limit for soft X-ray sources is 0.51 c s −1 at 0.3–2.0 keV, or 4.3 × 10 −13 erg cm −2 s −1 for a 100 eV blackbody. These limits would be four times lower if exploratory GTIs accumulate 10 ks of data after filtering at the level prescribed for faint sources. Such filtering selects background GTIs 85% of the time. An application of the model to a 1 s timescale makes it possible to distinguish source flares from possible surges in the background. 

Abstract Magnetars, isolated neutron stars with magnetic-field strengths typically ≳10 14 G, exhibit distinctive months-long outburst epochs during which strong evolution of soft X-ray pulse profiles, along with nonthermal magnetospheric emission components, is often observed. Using near-daily NICER observations of the magnetar SGR 1830-0645 during the first 37 days of a recent outburst decay, a pulse peak migration in phase is clearly observed, transforming the pulse shape from an initially triple-peaked to a single-peaked profile. Such peak merging has not been seen before for a magnetar. Our high-resolution phase-resolved spectroscopic analysis reveals no significant evolution of temperature despite the complex initial pulse shape, yet the inferred surface hot spots shrink during peak migration and outburst decay. We suggest two possible origins for this evolution. For internal heating of the surface, tectonic motion of the crust may be its underlying cause. The inferred speed of this crustal motion is ≲100 m day −1 , constraining the density of the driving region to ρ ∼ 10 10 g cm −3 , at a depth of ∼200 m. Alternatively, the hot spots could be heated by particle bombardment from a twisted magnetosphere possessing flux tubes or ropes, somewhat resembling solar coronal loops, that untwist and dissipate on the 30–40 day timescale. The peak migration may then be due to a combination of field-line footpoint motion (necessarily driven by crustal motion) and evolving surface radiation beaming. This novel data set paints a vivid picture of the dynamics associated with magnetar outbursts, yet it also highlights the need for a more generic theoretical picture where magnetosphere and crust are considered in tandem. 

Abstract We report on NICER X-ray monitoring of the magnetar SGR 1830−0645 covering 223 days following its 2020 October outburst, as well as Chandra and radio observations. We present the most accurate spin ephemerides of the source so far: ν = 0.096008680(2) Hz, ν ̇ = − 6.2 ( 1 ) × 10 − 14 Hz s −1 , and significant second and third frequency derivative terms indicative of nonnegligible timing noise. The phase-averaged 0.8–7 keV spectrum is well fit with a double-blackbody (BB) model throughout the campaign. The BB temperatures remain constant at 0.46 and 1.2 keV. The areas and flux of each component decreased by a factor of 6, initially through a steep decay trend lasting about 46 days, followed by a shallow long-term one. The pulse shape in the same energy range is initially complex, exhibiting three distinct peaks, yet with clear continuous evolution throughout the outburst toward a simpler, single-pulse shape. The rms pulsed fraction is high and increases from about 40% to 50%. We find no dependence of pulse shape or fraction on energy. These results suggest that multiple hot spots, possibly possessing temperature gradients, emerged at outburst onset and shrank as the outburst decayed. We detect 84 faint bursts with NICER, having a strong preference for occurring close to the surface emission pulse maximum—the first time this phenomenon is detected in such a large burst sample. This likely implies a very low altitude for the burst emission region and a triggering mechanism connected to the surface active zone. Finally, our radio observations at several epochs and multiple frequencies reveal no evidence of pulsed or burst-like radio emission. 

Analyses of pulsar timing data have provided evidence for a stochastic gravitational wave background in the nanohertz frequency band. The most plausible source of this background is the superposition of signals from millions of supermassive black hole binaries. The standard statistical techniques used to search for this background and assess its significance make several simplifying assumptions, namely (i) Gaussianity, (ii) isotropy, and most often, (iii) a power-law spectrum. However, a stochastic background from a finite collection of binaries does not exactly satisfy any of these assumptions. To understand the effect of these assumptions, we test standard analysis techniques on a large collection of realistic simulated data sets. The data-set length, observing schedule, and noise levels were chosen to emulate the NANOGrav 15 yr data set. Simulated signals from millions of binaries drawn from models based on the Illustris cosmological hydrodynamical simulation were added to the data. We find that the standard statistical methods perform remarkably well on these simulated data sets, even though their fundamental assumptions are not strictly met. They are able to achieve a confident detection of the background. However, even for a fixed set of astrophysical parameters, different realizations of the universe result in a large variance in the significance and recovered parameters of the background. We also find that the presence of loud individual binaries can bias the spectral recovery of the background if we do not account for them.

 

Abstract We report the discovery of the unusually bright long-duration gamma-ray burst (GRB), GRB 221009A, as observed by the Neil Gehrels Swift Observatory (Swift), Monitor of All-sky X-ray Image, and Neutron Star Interior Composition Explorer Mission. This energetic GRB was located relatively nearby ( z = 0.151), allowing for sustained observations of the afterglow. The large X-ray luminosity and low Galactic latitude ( b = 4.°3) make GRB 221009A a powerful probe of dust in the Milky Way. Using echo tomography, we map the line-of-sight dust distribution and find evidence for significant column densities at large distances (≳10 kpc). We present analysis of the light curves and spectra at X-ray and UV–optical wavelengths, and find that the X-ray afterglow of GRB 221009A is more than an order of magnitude brighter at T 0 + 4.5 ks than that from any previous GRB observed by Swift. In its rest frame, GRB 221009A is at the high end of the afterglow luminosity distribution, but not uniquely so. In a simulation of randomly generated bursts, only 1 in 10 4 long GRBs were as energetic as GRB 221009A; such a large E γ ,iso implies a narrow jet structure, but the afterglow light curve is inconsistent with simple top-hat jet models. Using the sample of Swift GRBs with redshifts, we estimate that GRBs as energetic and nearby as GRB 221009A occur at a rate of ≲1 per 1000 yr—making this a truly remarkable opportunity unlikely to be repeated in our lifetime. 

The radio galaxy 3C 66B has been hypothesized to host a supermassive black hole binary (SMBHB) at its center based on electromagnetic observations. Its apparent 1.05 yr period and low redshift (∼0.02) make it an interesting testbed to search for low-frequency gravitational waves (GWs) using pulsar timing array (PTA) experiments. This source has been subjected to multiple searches for continuous GWs from a circular SMBHB, resulting in progressively more stringent constraints on its GW amplitude and chirp mass. In this paper, we develop a pipeline for performing Bayesian targeted searches for eccentric SMBHBs in PTA data sets, and test its efficacy by applying it to simulated data sets with varying injected signal strengths. We also search for a realistic eccentric SMBHB source in 3C 66B using the NANOGrav 12.5 yr data set employing PTA signal models containing Earth term-only as well as Earth+pulsar term contributions using this pipeline. Due to limitations in our PTA signal model, we get meaningful results only when the initial eccentricitye₀< 0.5 and the symmetric mass ratioη> 0.1. We find no evidence for an eccentric SMBHB signal in our data, and therefore place 95% upper limits on the PTA signal amplitude of 88.1 ± 3.7 ns for the Earth term-only and 81.74 ± 0.86 ns for the Earth+pulsar term searches fore₀< 0.5 andη> 0.1. Similar 95% upper limits on the chirp mass are (1.98 ± 0.05) × 10⁹and (1.81 ± 0.01) × 10⁹M_☉. These upper limits, while less stringent than those calculated from a circular binary search in the NANOGrav 12.5 yr data set, are consistent with the SMBHB model of 3C 66B developed from electromagnetic observations.

 

The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has reported evidence for the presence of an isotropic nanohertz gravitational-wave background (GWB) in its 15 yr data set. However, if the GWB is produced by a population of inspiraling supermassive black hole binary (SMBHB) systems, then the background is predicted to be anisotropic, depending on the distribution of these systems in the local Universe and the statistical properties of the SMBHB population. In this work, we search for anisotropy in the GWB using multiple methods and bases to describe the distribution of the GWB power on the sky. We do not find significant evidence of anisotropy. By modeling the angular power distribution as a sum over spherical harmonics (where the coefficients are not bound to always generate positive power everywhere), we find that the Bayesian 95% upper limit on the level of dipole anisotropy is (C_l=1/C_l=0) < 27%. This is similar to the upper limit derived under the constraint of positive power everywhere, indicating that the dipole may be close to the data-informed regime. By contrast, the constraints on anisotropy at higher spherical-harmonic multipoles are strongly prior dominated. We also derive conservative estimates on the anisotropy expected from a random distribution of SMBHB systems using astrophysical simulations conditioned on the isotropic GWB inferred in the 15 yr data set and show that this data set has sufficient sensitivity to probe a large fraction of the predicted level of anisotropy. We end by highlighting the opportunities and challenges in searching for anisotropy in pulsar timing array data.

 

We have used X-ray data from the Neutron Star Interior Composition Explorer (NICER) to search for long-timescale temporal correlations (“red noise”) in the pulse times of arrival (TOAs) from the millisecond pulsars PSR J1824−2452A and PSR B1937+21. These data more closely track intrinsic noise because X-rays are unaffected by the radio-frequency-dependent propagation effects of the interstellar medium. Our search yields strong evidence (natural log Bayes factor of 9.634 ± 0.016) for red noise in PSR J1824−2452A, but the search is inconclusive for PSR B1937+21. In the interest of future X-ray missions, we devise and implement a method to simulate longer and higher-precision X-ray data sets to determine the timing baseline necessary to detect red noise. We find that the red noise in PSR B1937+21 can be reliably detected in a 5 yr mission with a TOA error of 2μs and an observing cadence of 20 observations per month compared to the 5μs TOA error and 11 observations per month that NICER currently achieves in PSR B1937+21. We investigate detecting red noise in PSR B1937+21 with other combinations of observing cadences and TOA errors. We also find that time-correlated red noise commensurate with an injected stochastic gravitational-wave background having an amplitude ofA_GWB= 2 × 10⁻¹⁵and spectral index of timing residuals ofγ_GWB= 13/3 can be detected in a pulsar with similar TOA precision to PSR B1937+21. This is with no additional red noise in a 10 yr mission that observes the pulsar 15 times per month and has an average TOA error of 1μs.